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Abstract 

This paper is a contribution to the theory of coherent crystals. We present arguments claiming that 
negative minima in the Fourier transform of a soft pair interaction may give rise to the coexistence 
of diagonal and off-diagonal long-range order at high densities, and that this coexistence may be 
detectable due to a periodicity seen on the off-diagonal part of the one-body reduced density 
matrix, without breaking translation invariance. As an illustration, we study the ground state 
of a homogenous system of bosons in continuous space, from the interaction retaining only the 
Fourier modes v(k) belonging to a single nonzero wave number |k| = q. The result is a mean-field 
model. We prove that for w(k) > the ground state is asymptotically fully Bose-condensed, while 
for w(k) < at densities exceeding a multiple of h 2 q 2 /2m|w(k)| it exhibits both Bose-Einstein 
condensation and diagonal long-range order, and the latter can be seen on both the one- and the 
two-body density matrix. 

PACS: 03.75.Hh, 05.30.Jp, 67.80.K- 



1 Introduction 

Ever since the clarification by Penrose and Onsager [1] of the notion of Bose-Einstein condensation 
(BEC) for interacting bosons, an important question has been the accommodation of a diagonal long- 
range order (DLRO) with the off-diagonal one (ODLRO), characterizing BEC. This question has 
gained a new interest with the experimental discovery of nonclassical rotational inertia in solid helium 
[2] , its being a possible signature of coexistence of the two types of order [3] . Penrose and Onsager 
defined BEC as the non-decay of Ajv /N with N going to infinity, where N is the number of particles and 
A at is the largest eigenvalue of the one-body reduced density matrix. If the Hamiltonian is translation 
invariant, as in the case when it contains only a shift-invariant interaction and no external field and 
is defined with periodic boundary conditions on a cube A, a Perron-Frobenius argument shows that 
Xn = (o-Qao), the expectation value of the occupation number of the zero-momentum one-particle 

state. Expressing Oq*^ in terms of the field operators aM(x), one obtains Ajy = / A (a*(x)a(0)) dx. 
This, together with the positivity of (a*(x)a(0)), immediately yields that BEC is equivalent to the 
non-decay of (a*(x)a(0)) as |x| — ■> oo (after the thermodynamic limit), that is, ODLRO - a conclusion 
obtained already by Penrose and Onsager [4]. 

In contrast to ODLRO which has no classical analogue, DLRO refers to classical spatial ordering 
of particles. The most common form of this is the periodic order of a crystal. Being classical, DLRO 
should be better understood than ODLRO, but this is not the case. Because of the enormous difficulty 
to deal with structural ordering in continuous space, we have only a few very specific rigorously proven 
examples of phase transitions [5] - [7] and of ground-state ordering ([8] - [10] are examples with relevance 
for the forthcoming discussion), all about classical systems and most of them relatively recent. A part 
of the new activity [9] - [11] concerns the classical analogue of the so-called coherent crystals, a special 
kind of supersolids, which is the subject of the present paper. The name and a detailed study appears 
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in Refs. [12, 13], but the origin is earlier, in Vlasov's theory of solids based on the Boltzmann equation 
[14] and in subsequent descriptions of the ground state of fermion or boson systems in a self-consistent 
single-particle (Hartree or Hartree-Fock) approximation [15] - [17]. While Vlasov's coherent crystal can 
be classical, that of the other papers is an intrinsically quantal object resulting from the crystallization 
of a quantum fluid. The conclusion of these studies is that for certain interactions at high density 
the expected single-particle ground state is a periodic function, therefore the particle density exhibits 
spatial periodicity and, hence, DLRO. For the authors of Ref. [12] a coherent Bose crystal necessarily 
displays BEC, so it is a supersolid. If we disregard the approximations leading to it, the defining 
characteristic of a coherent crystal is that its lattice constant is not determined by the particle density 
but by the global properties of the interaction. As a result, when the density is increased, the lattice 
constant is unchanged and the number of particles within the unit cell increases. Moreover, in the 
quantum crystal the kinetic energy per particle is considerably larger than the interaction energy per 
pair [13], which can hold only if the particle mass is small and the interaction is weak. Thus, the 
particles move relatively freely but their motion must be coherent, otherwise they could not create 
and maintain an effective periodic field. The existence of coherent quantum crystals has remained 
hypothetical, but their classical counterpart has recently been recognized to exist in the soft matter 
world [11]. Below the term "coherent crystal" will be applied to both the classical and the quantal 
object. The coherent quantum crystal is understood to exhibit both DLRO and BEC. "Coherent 
supersolid" will be used as its synonym. 

One can fully appreciate the peculiarities of coherent crystals only by a comparison with "ordinary" 
crystals composed of identical particles. Ordinary crystals are thought to consist of localized particles 
making oscillations about separate lattice sites, their respective equilibrium position. The properties 
of ordinary crystals can be well described by a model based on this image, assuming harmonic (in 
the Born-Karman theory) or anharmonic restoring forces [18, 19]. Thus, each particle is tagged by a 
lattice site, and the pair potential docs not depend on the difference of the position vectors but on 
the difference of the displacements from equilibrium. This means a built-in violation of permutation 
invariance. At a certain strength of anharmonicity the restoring forces vanish, the particles get freed 
and the model breaks down. If an analytic continuation beyond the singularity were possible, one 
could probably see the emergence of a coherent crystal. 

We may say that the principal difference between coherent and ordinary crystals is revealed by 
their fundamentally different approximate theories outlined above. To see also the similarities, a more 
adequate approach would be to work in both cases with a full many-body Hamiltonian containing 
a permutation-invariant interaction. Then, one would never see the spontaneous breakdown of the 
permutational symmetry in an equilibrium state (apart from the cases of classical ground states, 
compressed hard-core systems or systems with unstable interactions [6]). Only translation invariance 
can be broken, and the distinction becomes subtle, with a gradual transition between ordinary and 
coherent crystals as the interaction varies. Also, a system can show the properties of a coherent crystal 
at lower densities and those of an ordinary crystal at high densities. The prototype of an ordinary 
solid is presumably obtained with a purely repulsive interaction having a nonintegrable divergence 
at the origin. One then expects that the formation of a crystal at high density will be governed by 
short-range correlations, the lattice constant will be determined by the density and the crystal will 
contain a single particle in its unit cell. The prototype of a coherent crystal is provided by bounded 
interactions with a partly negative Fourier transform. Classically this is now confirmed [11], and the 
works [12, 13] are based on such an interaction. The repulsive square-well interaction, mentioned 
explicitly by Vlasov [cf. his equation (55)], has this property and defines a clear example of a classical 
coherent crystal [20]. 

An important intermediate case is an interaction with a strong short-range repulsion and an 
attractive part at a characteristic distance. Such interactions appear in Refs. [14, 16, 17]. The best 
known example is the Lennard-Jones potential, introduced to describe the Van der Waals interaction 
between neutral atoms and, therefore, relevant in the treatment of crystals of noble gas atoms. For He 4 
and He 3 the zero-point kinetic energy can be large compared with the interaction energy of a pair, cf. 
[17], but the decoupling of the lattice constant from the density may be satisfied only approximately. 
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Because of the negative minimum in the pair potential, the interaction prefers a given structure of a 
well-defined density p s , in the sense that the classical ground state energy per volume as a function 
of the particle density p attains its absolute minimum on this structure. If p ^ p s but is close to it, 
\p — p s \ <C p s , without a minimum at p s one would expect a crystal of density p. For the interaction in 
question, however, the lattice density may not adjust to that of the particles but follow it with some 
delay, e.g., it may be p s + c(p — p s ) 3 where < c < (p — p s )~ 2 - If this indeed occurs then the classical 
ground state, although still showing DLRO, is not completely ordered, there are particles or holes 
in excess with respect to a perfect crystal, leading possibly to a macroscopic zero point entropy. In 
the quantum mechanical ground state these particles or holes may condense and can be the source of 
ODLRO accompanying DLRO [17]. One must stress, however, that the actual knowledge about the 
classical ground state is insufficient to confirm the above scenario. For a two-dimensional analogue 
of the Lennard-Jones potential there is a proof of crystallization in the preferred structure [8] (the 
lattice constant is close to the distance at which the pair potential reaches its minimum) , but nothing 
is known rigorously for p ^= p s . 

Some of the earliest speculations about the coexistence of ODLRO and DLRO were based less on 
the properties of the Hamiltonian than on the supposed form of the ground state. Taking into account 
the permutation (but not the translation) invariance of the Hamiltonian, the quantum mechanical 
ground state of a perfect ordinary crystal of identical bosons would be 

N 

cD(ri,..., rjv ) =C^I]^( r -W- R 0- (!) 

Here C stands for normalization, the summation is over the permutations 7r of 1, . . . , N, the product 
goes over the vertices Ri of a lattice in A, and is a real nonnegative L 2 -function having a maximum 
at the origin and being practically independent of A for large volumes. The existence of such a ground 
state for a realistic Hamiltonian is as hypothetical as the existence of an ideal coherent quantum 
crystal whose ground state would be 

N 

*(r 1 ,...,r JV ) = JJV(ri). (2) 

i=l 

Here ip is a real nonnegative, nonconstant periodic function. It is easy to check that in the state $ 
there is no ODLRO while in ^ there is. Penrose and Onsagcr assumed that a spatially ordered ground 
state can only be a perfect crystal, and naturally concluded that there can be no BEC in such a state. 
On the other hand, the Hartree approximation in [12] naturally lead to a coherent supcrsolid state of 
the form (2). A more subtle A^-body wave function with probable ODLRO and DLRO was introduced 
by Chester [21], whose work prompted Leggett's suggestion [3] about the existence of a nonclassical 
rotational inertia in Bose solids. Apart from [12], studies of the coexistence of ODLRO and DLRO 
based on the properties of the Hamiltonian are mainly restricted to lattice models. Supersolid phases, 
in the sense of the coexistence, were found on bipartite lattices away from half-filling by approximate 
analytical [22, 23] and numerical [24] - [26] methods. 

This paper is meant to be a contribution to the theory of coherent quantum crystals. Even though 
their study has been neglected in the last years in favor of lattice models, they can play an important 
role in the explanation of the properties of light-atom quantum crystals at moderate pressures when 
the average interaction energy per pair is small compared with the kinetic energy per particle. An 
additional motivation comes from the recent discovery of classical coherent crystals in soft matter 
physics [11]. Based on earlier works [12, 13, 11], the appropriate interaction for the description of 
coherent crystals is bounded and has a partly negative Fourier transform. Without being able to deal 
with a full interaction, we go beyond the Hartree approximation [12, 13] by obtaining results for a 
mean-field interaction which contains only Fourier components of a given wave length. 

In Section 2 we recall the results on classical ground states for soft interactions and propose a 
method to detect the coexistence of ODLRO and DLRO in coherent quantum crystals. Section 3 
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contains the proof that if the Fourier transform of the interaction is partly negative, in high-density 
ground states of Bose systems the static structure factor develops a maximum on the support of the 
negative part. It is in Section 4 that we introduce and discuss the mean-field model. We prove that 
for a positive Fourier transform the ground state is fully Bosc-condcnsed at all densities while if the 
interaction has a negative Fourier transform then at high enough densities the ground state displays 
joint DLRO and ODLRO. The paper ends with a few remarks in Section 5. 

2 A signature of coexisting diagonal and off-diagonal long- 
range orders 

In relation with the Bose gas bounded or more generally integrable pair interactions are interesting 
as being at the heart of Bogoliubov's theory. Indeed, this theory is based on the Hamiltonian 



where A = [0,L] d , a d-dimensional cube of volume V = L d on which Ha is defined with periodic 
boundary conditions, A* = is the dual of A, e(k) = ti 2 k 2 /2m, AT k = a k a k , 



and a k , a k are the usual boson creation and annihilation operators. The underlying hypothesis in (3) 
is that v(k) is the Fourier transform of an integrable pair potential u(r) and is, therefore, a bounded 
continuous function. An integrable u, even if unbounded, is quantum mechanically a soft interaction 
in the sense that the potential energy in a two-particle state 4 r (ri, r2) may be finite even if ^(r, r) ^ 0. 
If v is also integrable then u is bounded and, hence, soft also classically. 

Ordering of classical soft particles exhibits some peculiar features. This can already be seen in 
the so far unique case treated rigorously [9, 10] when v(k) > and identically vanishes above some 
wave number. Above a threshold density the classical ground state is highly degenerate, the energy- 
minimizing configurations are continuously deformable and their unions (even overlapping ones) are 
also energy-minimizing at a higher density. This is the manifestation of the striking property that any 
ground state configuration creates a force-free field on a test particle. According to a recent result 
by Likos et al. [11], stacking of particles in classical ground states should be even more pronounced 
if v(k) has a negative minimum: this minimum determines the lattice constant of a crystal, and the 
particles occupy the lattice points independently of the density [27]. This seems to indicate that the 
crystal creates a periodic field on a test particle with minima in lattice points, even though the pair 
interaction may be purely repulsive. What these authors found is, in effect, the classical coherent 
crystal realized in soft matter physics. It is to be noted that the results of Ref. [11] were obtained 
by classical density-functional theory, an analogue of the Hartree approximation used in [12]. A hint 
about the role of a negative minimum in v(k) is provided by the Poisson summation formula 



written in a symmetric form. Here B is a Bravais lattice of density ps and B* is its reciprocal lattice. 
On the basis of this formula it seems reasonable that for minimizing the energy density the lattice 
constant of B or B* should be adjusted to the minimum of u or v, respectively. In the second case 
the selection of B is indirect, the direct selection occurring in dual space via the dual interaction v. 

In view of the above result on classical coherent crystals the conjecture that the ground state of 
the Hamiltonian (3) with a partly negative v is a coherent supersolid, showing both BEC and periodic 
order, seems natural. Below we argue that the coexistence of the two kinds of order can be detected 




(3) 




(4) 



k,k'GA* 




(5) 
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solely on the off-diagonal part of the one-body reduced density matrix, without breaking translation 
invariance. Concerning the periodic order the suggestion is new. To mark this fact, we shall use 
the term SLRO (spatial or structural long-range order) instead of DLRO which can be seen on the 
diagonal part of the two-body density matrix. Let us recall the relevant definitions. In a translation- 
invariant setting the iV-particle density matrix /5a, n (= |3>a)(3>a| at zero temperature, where <1>a is the 
TV-particle ground state of Ha) commutes with the total momentum operator and so do its reduced 
density matrices. With the short-hand (A) = Tr /5a, at A we then obtain 

(a*(x)a(y)) = (a*(x - y)a(O)) = £ (A^e"^-^ (6) 

keA* 

for the one-body reduced density matrix, and 

(a*(x)o*(y)a(x)a(y)) = V~ 2 £ (SOe"*^) (7) 

keA* 

for the diagonal part of the two-body reduced density matrix. Note that (£k) = iV[5(k) — 1] where 5(k) 
is the static structure factor. To detect long-range order, in the above formulas the thermodynamic 
limit must be taken. Below, lim without subscript will mean V, N -> oo and N/V = p. Let o-W(x-y) 
and cr( 2 )(x — y) be the respective limits, admitting they exist. Both are real nonnegative functions 
[28], the Fourier transform of is also positive so that a^(x) < a^\0) = p, and the average of 
<T^(x) is p 2 . Now DLRO is seen on the non- decay of (x) - p 2 and ODLRO is seen on a^ 1 ' (x) -t* 
and on its non- vanishing average, 



v = lim (2s)~ d / cr (1) (x) dxi . . . dx d > 0, 

s ^°° J\xi\<s 



(8) 



giving the condensate density. 

In the thermodynamic limit (N^) /V tends to some finite positive measure dv(k) on M d . The 
general form of crW(x) is therefore 

a«(x)= J e- ik - x dKk)=a«(x)+4 1 p )(x)+a( c 1 )(x), (9) 

the three terms corresponding to the decomposition of d^(k) into absolutely continuous (ac), pure 
point (pp) and singularly continuous (sc) parts. For (9) to hold, J di/(k) = p must be satisfied, thus, 

F(ko) = I diz(k) = lim J2 { -ir^p ( 10 ) 

• 7 l k l <fco keA*,|k|<fe 

must tend to p as k — > oo. This is indeed the case, otherwise the kinetic energy density would diverge 
in the thermodynamic limit because of the inequality 

t = lim £ e(|k|)^ > e(fco) [p F(k )} . (11) 

keA* 

Actually, F(ko) > p — c/k 2 with some c > as fco — > oo. 

Due to the Riemann-Lebesgue lemma, a^(x) tends to zero as |x| goes to infinity. This is not true 
for Opp-'(x) and uic^x). The pp and sc parts of dv are concentrated on sets of zero Lebesgue measure 
and can be interpreted as different kinds of (generalized) BEC [29]. However, the pure point part 

4 1 p )(x)=z, + Yl ^e- ik - x ( 12 ) 
kezX 1 ' 

where is a countable subset of W 1 not containing has another, more relevant interpretation. In 
[30] it was shown that if > for some k ^ then there exist infinite-volume equilibrium states 
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which are periodic in the direction k. More generally, the sum in Eq. (12) can describe almost periodic 
order if is a general nonempty set, and £>-periodic order if = B* . Whether or not ail\x) 
can be nonzero is unknown, but in principle it can also contribute to the asymptotic (|x| — > oo) 
non-constancy of crW(x). In the definition below we disregard this possibility. 

Definition. We say that a density matrix displays spatial or structural long-range order (SLRO) 
if its associated set defined in Eq. (12) is nonempty. 

From the definition it follows that SLRO implies BEC. Indeed, if v-^ > for some k ^ then 
vq > as well: 

Vk = lim (2s)- d { e ik x cr (1) (x) Ax x . . . dx d < v (13) 

because > 0. Furthermore, SLRO implies also DLRO. This is an interpretation of the result 
of Ref. [30] and follows from the meaning of = t/_k > 0: If translation invariance were broken 
e.g. by an arbitrarily weak k-periodic external field, one would see the macroscopic occupation of a 
one-particle state ~ cos(k • r + a). This would imply the spatial periodicity of the density and would 
show up in the diagonal part of the one-body density matrix. Hence, we conclude that SLRO implies 
the coexistence of ODLRO and DLRO. 



3 Effect of a partly negative v(k) on the ground state of a 
system of bosons 

In this short section we show that whenever the Fourier transform of the interaction has a negative 
part V-, at high enough densities the static structure factor 5(k) develops a maximum on the support 
of V-. Decompose v into positive and negative parts, 

u(k) = v+(k)-v_(k) (14) 

where v± > and v + v- = 0, and suppose that V- is not identically vanishing. Because u € L 1 (K <i ), v 
is a bounded continuous function decaying at infinity, therefore v(k) < on a union of open intervals. 
Choose a k such that v(k ) < 0, and compute the energy in the product state 

N [2 1 

*(r)w=nvy cos 2 k °' rj - (15) 

Here (r) N = (n, . . .,r N ). For 



H A = H A -^N(N-1) (16) 



by a simple computation one finds 



i(^ A |*) = i[e(|ko|)+Mko)]p. (17) 
This is an upper bound to the ground state energy density, 

^^ A \H A \^ A )=t + ^ £ «(k)[5(k)-l]<i[6(|ko|) + P«(ko)]p. (18) 

0#keA* 

Insert here (14), use the lower bound for 5(k) on suppu+ and also for t to find 

E keA . ^-(k)[S(k) - 1] P Kko)| - 6(kp) - 1 E 0#keA . Mk) 

E k6A .«-(k) - fE k£A .»-(k) ■ 1 ' 

Note that v(0) — Ju(r)dr > and therefore v(k) > in a neighborhood of k = 0. The weighted 
mean of S(k) — 1 on the left-hand side of (19) becomes positive if p is above a volume-independent 
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threshold value and then it increases at least linearly with p. Since the arithmetic mean of S(k) over 
A* \ {0} is 1, wc can indeed see the appearance of a maximum of 5(k) somewhere on suppw_. This 
still docs not prove DLRO — for that we should show that 5(k) <~ N for some k ^ at a fixed 
density p as V, N — > oo , N/V — p — but indicates its possible occurrence. 

4 A mean-field model for coherent supersolids 

Unfortunately, the present mathematical knowledge does not allow one to prove even BEC alone for 
the Hamiltonian (3) with a general v(k). However, we can check the ideas outlined in Section 2, in 
particular the role of a negative minimum of v(k) in a simplified case. Let us recall that 

v(k) = / u(r)e- 4kr dr= /'w A (r)e- lkr dr, (20) 

JM d J A 

where 

u A (r) = U ( r + Ln ) ( 21 ) 
nez d 

is the periodized interaction that we normally use when we work with periodic boundary conditions. 
By inverting (20), 

%(r) = \ £ v(k)e^. (22) 

keA* 

The constant first term, v(0)/V, is known as the mean-field interaction for Bose gases. Although in 
the canonical ensemble the corresponding mean-field model does not differ from the ideal Bose gas, 
the mean-field interaction is superstable (the potential energy of N particles is v(0)N(N — 1)/2V), 
which makes it possible to treat the system in the grand-canonical ensemble for arbitrarily large values 
of the chemical potential and, hence, in the region of BEC. The mean-field model has been studied 
extensively in the past; rigorous results can be found e.g. in Refs. [31]-[33]. The mean-field interaction 
is interesting also as the limiting case of a scaled interaction (the Kac potential); this observation was 
the basis of recent attempts [34, 35] to prove BEC for nontrivial interactions. If we go further in the 
expansion (22) by retaining a finite number of terms, the result is still a mean-field interaction which 
is, however, not constant. By a proper choice of the additional terms one may hope to preserve at least 
some qualitative features of the fully interacting system. To demonstrate that a negative minimum 
in v (k) may lead to joint DLRO and ODLRO the proper choice is the set of wave vectors of length q 
at which the minimum is attained. For a cubic domain this leads to the interaction 

. . a 2b , . 

«a(x) = — + — 2^ cosc b ( 23 ) 

where = qej. Let L be an integer multiple of 2n/q, then A = [0, L] d is a period cell for the Bravais 
lattice B = (2ir/q)Z d , and A* = {2ir/L)Z d D B* = q7L d . The simple cubic B is just an example. 
With suitable q 3 and A, B can be any Bravais lattice. If a > 2d\b\ then ma is superstable. However, 
a does not affect the eigenstates of H\, and the potential energy density at a fixed particle density 
is bounded below without the term a/V. We can therefore set a = 0, and for comparison we shall 
consider both b > and b < 0. As the theorem below will show, for b < (23) defines the mean-field 
model of coherent supersolids in continuous space. 

Before proceeding further, let us comment on the possible relevance of the mean-field model in the 
study of the Hamiltonian (3). A mean-field-type truncation (23) of a soft interaction is a high-density 
approximation. If v > 0, there is a clear indication for this being true. In effect, from Eqs. (16) and 
(18), and by using fi(r)jv = V~ N / 2 as a variational wave function instead of (15), one obtains the 
bounds 

v(0) 2 u(0) ,. 1 _ . TT _ , v(0) 2 

^-V - ~y-P < e g ., - lim-(<D A |iJ A |$ A ) < A2p 2 (24) 
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for the ground state energy density (see also the Appendix of Lieb [36]). Thus, as p tends to infinity, 
the leading term in e g . s . is v(0)p 2 /2, the exact result for the mean-field Bose gas. If v is partly 
negative, the precise asymptotic (p — > oo) form of the ground state energy density is unknown. It is 
clear, however, that minimization of the potential energy dictates that if the classical ground state is 
ordered, this order shows up in the quantum mechanical ground state. Our discussion in the previous 
section supports this expectation. The result obtained by Likos et al. [11, 27] reveals the order 
characteristic to classical coherent crystals, particles accumulating on the sites of a given lattice in 
high-density ground states. It will be seen that the mean-field interaction (23) has such ground states 
for b < 0. Let us note that in Ref. [12] explicit mention is made of a pair potential of the form (23) as 
an interaction that in Hartree approximation would lead to a ground state in which the only possible 
momentum transfer is for |k| = q. For the full many-body ground state this view will prove to be 
over-simplified. Finally, we stress that although with (23) a periodicity is introduced by hand, the 
interaction is still translation- invariant, a case not to be confounded with the presence of a periodic 
external field as in optical lattices [37]. Moreover, the periodicity will be seen to disappear in the 
thermodynamic limit if b > and to stay if b < and the density is high enough. The difference 
between the two cases is the main point of the following theorem. 

Theorem. For b > the ground state of H A is asymptotically fully Bose-condensed at any p. For 
b < and p > (e/4)e(g)/|6| the ground state presents both BEC and SLRO/DLRO in the sense that 
both and are nonconstant i?-periodic functions. 

Proof. 1. We start with identifying the classical ground state configurations and energies. The 
TV-particle interaction energy is 

Ua(t)n = M r * - r j) = ^ 

For b > 0, 

U A (r) N > -dbN/V (26) 

and the lower bound is actually the classical ground state energy. Any (r)jv making J2n=i e * qj " r " 
disappear for all j is a classical ground state. Each of these (r) jv can be obtained by the following 
construction. Writing r„ = (2-k /q)(x\ n , . . . ,Xd n ), • r„ = 2nXj n . For each j we partition the 
coordinates {xj n }n=i m t° groups of at least 2 elements and in a group of m elements we choose a real 
a and assign the values a, a+l/m, . . . , a+(m— l)/m to the coordinates in an arbitrary order. The sum 
of the exponentials will disappear separately for each group. In this highly degenerate set of ground 
state configurations only a negligible number shows long-range order. All will be mixed in a unique 
unordered quantum mechanical ground state; that this one is asymptotically fully Bose-condensed, is 
not a priori obvious. 
For b < 0, 

U A (r) N > ~d\b\N 2 /V + d\b\N/V (27) 

and the lower bound is again the classical ground state energy. All (r)jv giving e lqi ' r ™ = 1 for each 
j and n and all translates of such configurations are classical ground states. Independently of N, 
these are the (r)jv with r„ — r m G B for every m,n. Thus, for b < in any classical ground state 
configuration all particles occupy the points of a lattice whose lattice constant is independent of the 
particle density. Classically, the simple interaction (23) reproduces the zero-temperature result of 
Likos et al. [11, 27]. 

2. The Hamiltonian (3) with 

d 

v(k) - a<5 k , + &X(<5k, qj + <5k,- qj ) (28) 
i=i 

corresponding to (23) has a large number of extra symmetries. Consider the coset group A* / B* of 
A*. For any C G A* / B* , H A commutes with Nc — X)kec^ k - Consider now any partition of A, 



N 



71=1 



dbN 



(25) 
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1 < ni < ■ ■ ■ < n p < N, J2 n i — ^ an d any C\, . . . ,C P e A* / B* . Then the vectors * satisfying 
— for j = l,...,p form an iJA-invariant subspace of [L 2 (A) N ] sym . For the iV-particle 
ground state $a of H\, Nb*$a = N&\, meaning that $a is constructed solely from plane waves 
whose k vector is in B*. Indeed, <&a is the eigenstate of e~@ HA belonging to the largest iV-particle 
eigenvalue. (Here (3 is an arbitrary positive number of dimension 1/energy.) Because the integral 
kernel of e~@ HA is positive, $a is positive and hence unique. Therefore, it must show all symmetries 
of H\. Thus, <I>a is permutation and translation invariant, and B-periodic separately in each position 
vector. The only invariant subspace whose elements have these symmetries is Nb*^ — 

3. In the ground state, (•) = ($a| • |$a) and (6) reads 

(a*(x)a(y)) =N f $ A (x, (r)f ) $ A (y, (r)?) d(r)f . (29) 

According to point 2, 

(a*(x)a(y)) = V' 1 ]T (N^e-^^. (30) 

So (a*(x)a(0)) is 5-periodic in finite volumes. By the diagonal process we can choose a subsequence of 
L values such that (iVk) /V converges to some for every k e B* as V and N = [pV\ tend to infinity. 
We have then three possibilities, (i) a^\x) = p = Vq, meaning full BEC. (ii) X)kes* = p > vq. 
Because < i/ , there is partial BEC (z^o > 0) and partial SLRO (fk > for some nonzero keB*), 
realized as a sort of full generalized BEC, a^{x) = <7 P p (x) is a nonconstant £?-periodic function, (iii) 
< X)kes* < Pi which can immediately be excluded because it would lead to a diverging kinetic 
(and total ground state) energy density, cf. Eq. (11). In cases (i) and (ii) the infinite sum and the 
limit are interchangeable in hm^ keB » (Nk)/V(— p). 

4. Let b > 0. Using Eqs. (25), (26) and fi(r)jv = V~ N I 2 as a variational wave function, we obtain 

-dbN/V = mmU A (r) N < ($ A |(7 A |$ A ) < ($a|#a|$a) = E g . B . < {Sl\H A \il) = 0. (31) 

(r)w 

The last equality holds because both the kinetic and the potential energies vanish in the state fi. 
Dividing by V and taking the thermodynamic limit this yields e g . s . = lim($A|-ffA| < l>A)/U = 0, and the 
kinetic and potential energy densities separately vanish in the ground state. On the other hand, 

t= e(|k|H > e(<z) ^ = <D(P-V°)- (32) 

kGB* O^keB* 

If vq < p, we would obtain t > 0; thus, v Q = p which is case (i) of point 3. This proves the b > part 
of the theorem. 

5. Let b < 0. First we show that is a nonconstant periodic function if the density is high 
enough. The ground state energy density reads 

e g .s. =t- \b\dps q (33) 

where 

8 q = VmS{qj)/V (for any j). (34) 

With the variational ansatz fi and Eq. (27) we find — d\b\p 2 < e g . s . < 0. A better upper bound is 
obtained with the trial function 

N 

*(r) Ar =n^(rn), VW = (9 2 /27r) d/4 e-« V2 / 4 (35) 

n=l 

where the parameter of the Gaussian is already optimized. This yields 

e g . s .<dp[e(q)/4-\b\p/e] (36) 
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which is negative if p > (e/A)e{q)/\b\. Thus, s q > for p > (e/4)e(g r )/|6|, and the expansion (7) (which 
also reduces to B*) reads 

d 

<jW (x) = p 2 + 2ps q cos • x + • • • (37) 
i=i 

To prove the periodicity of we start with Eq. (4), apply the Schwarz inequality to the addends, 
the Cauchy inequality to the sum and the equality X^keA* = ^ x ^he resu ^ is 

S(k)<l + N-(^(Noj-^(N^y , (38) 

holding in all generality. Dividing by V and letting N and V tend to infinity, for k = at p > 
(e/4)e(q)/\b\ we find 

(V^o-V 1 ^) <P- s q<P- (39) 
Now either v Q > > or = and then v < p follows from (39). Because X^keB* = P, 
v o > > for some nonzero keB*. This ends the proof of the theorem. 



5 Concluding remarks 

1. We have not proven that for b < there is no DLRO at low densities, but based on Eq. (36) 
the existence of a p c ~ e(q)/\b\ seems plausible. The order of the quantum phase transition at p c is 
an interesting open question. In Refs. [12, 13] there is a discussion of this point. It is argued that 
whenever there exists a triplet of reciprocal lattice vectors such that |kj| = q and ki + k2 +k3 = then 
by Landau's argument coherent crystallization is a first order transition. This would be the case if, 
instead of a simple cubic B, we would choose B to be the basic-centered cubic lattice or the hexagonal 
or rhombohedral lattice with some specific aspect ratios. On the other hand, if no such triplet exists 
- this is the case e.g. of a simple cubic B - then, because the interaction (23) is anisotropic, the 
transition could be of second order, cf. [13]. 

2. It may seem curious that DLRO is proven in a translation invariant situation, therefore the state 
in infinite volume is not a pure one but a uniform mixture over translations. However, the same holds 
true for the proof of BEC: by not breaking gauge invariance, the infinite- volume state is a uniform 
mixture over the phases of condensates [38] - [40] . 

3. In our proof we have not used that SLRO implies DLRO but proved the periodicity of and 
separately. It should be possible to prove the above implication in the general case without referring 
to the breakdown of translation invariance, by showing that the asymptotic periodicity of implies 
that of a< 2 ). 

4. The way $A(r)iv describes 5-periodic order is somewhat special, since it is -B-periodic separately 
in each r„. In general for a translation-invariant ground state one can expect only some weaker 
property, e.g. it having maxima at configurations (t)n such that r„ — r m e B for a positive fraction, 
nonvanishing in the limit of infinite volume, of pairs of particles (and r m = r„ excluded for interactions 
that are non-integrable at the origin). 

5. The possible appearance of SLRO in the one-body reduced density matrix for a general soft 
interaction is better seen on the ground state formula (29). Suppose that $a describes i3-periodic 
order in the weak sense given above. If x — y € B, (r)^ can be chosen in many different ways so 
that both $a (x, (r)/y) and <1>a (y, (r)^) are at a maximum. Then the integral itself may be at a 
maximum. 

6. Although the Poisson duality (5) places u and v on equal footing, quantum mechanically this 
symmetry is broken by the kinetic energy term. This is the reason why we expect in general different 
results when u has negative minima but v is positive and when v has negative minima and u > 0, 
even though in both cases there can be a lattice of a specific density minimizing the potential energy. 

7. The mean-field result suggests that in the ground state of the Hamiltonian (3) there is no DLRO 
without BEC. This may not be the case at positive temperatures. A crystal may be formed at a 
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temperature higher than the critical temperature of BEC at the given density. If now the temperature 
is kept fixed and the density is increased, the systems responds as a classical coherent crystal, by the 
accumulation of particles in each unit cell. However, increasing the density at a given temperature 
leads also to BEC: in this respect a system with a bounded interaction should not differ from the 
ideal Bose gas. Passing the critical density for BEC marks the transition from classical to quantum 
coherent crystal. 

Acknowledgement. I thank the referee for drawing my attention to some publications, in particular 
to references [12] and [13]. This work was supported by OTKA grants T 46129 and 67980K. 
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